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In  this  research  the  problem  cF  designing  a  controller 
for  the  dive  maneuver  of  ar.  Autonomous  Underwater  Vehicle 
(AUV)  is  addressed.  The  highly  nonlinear  nature  of  the  vehicle 
dynamics  and  the  requirement  for  fast  maneuvering  call  for 
robust  control  techniques.  In  particular,  Variable  Structure 
Control  (VSC)  combined  with  Adaptive  Control  (AC)  techniques 
seem  to  yield  satisfactory  performance  in  terms  of  robustness, 
capability  to  adjust  to  different  operating  conditions,  and 
speed  of  response.  Also,  linear  robust  techniques  based  on  LQG 
and  robust  observers  are  presented  to  address  the  case  when 
the  whole  state  (in  terms  of  pitch  rate,  pitch  and  depth)  is 
not  available  for  measurement. 
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I.  INTRODUCTION 


In  the  last  few  years,  considerable  interest  in  Autonomous 
Underwater  Vehicles  (A’JV)  has  arisen.  Possible  utilization 
scenarios  for  an  untethered  submersible  have  also  grown  and 
include:  sea  floor  mapping,  target  identification,  remote 
reconnaissance,  object  recovery,  etc.  A  typical  mission  would 
include  downloading  of  the  mission  program,  remote  positioning 
and  launch,  subsequent  vehicle  recovery,  and  data  collection. 
AUV  vehicles  now  in  design  or  production  can  employ  basically 
any  shape  from  a  box  shape  suited  to  offshore  work  to  a  body 
of  revolution  suited  for  high  speed  maneuverability  or  long 
distance  missions.  The  AUV  group  at  the  Naval  Postgraduate 
School  (NPS)  has  chosen  the  shape  shown  in  Figure  1  which  was 
derived  from  the  Swimmer  Delivery  Vehicle  (SDV)  [Ref.  1].  Part 
of  the  rationale  was  that,  since  extensive  data  exists  on  that 
particular  design,  it  would  facilitate  calculation  of  a 
reasonably  accurate  analytical  model  for  a  similar  but  smaller 
scale  model.  This  is  particularly  advantageous  to  the  design 
of  a  control  system  for  the  7*UV,  since  knowledge  of  an 
accurate  dynamic  model  is  in  general  the  basis  of  a  reliable 
control  design.  Also,  tests  conducted  by  computer 
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simulations  yield  results  which  are  closer  to  the  physical 
performance  of  the  vehicle. 

Ths  initial  goal  of  the  AUV  group  is  to  build  a  working 
prototype  AUV  that  will  be  able  to  expand  its  capabilities  to 
meet  the  increased  complexity  of  future  missions. 


As  with  all  analytical  models  of  submersibles,  the 
equations  of  motion  are  highly  nonlinear  and  vary  with  speed 
and  the  amount  of  equipment  onboard.  Adding  to  that  is  the 
difficulty  in  obtaining  exact  dynamic  coefficients.  Because 
of  the  high  level  of  nonlinearities  and  model  uncertainties, 
we  have  to  consider  more  sophisticated  design  techniques. 
Recently,  an  adaptive  controller  has  been  designed  by  Schwartz 
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[Ref.  2]  where  the  parameters  of  a  compensator  are  adjusted 
on-line  on  the  basis  of  an  estimated  transfer  function.  This 
approach  has  not  only  provided  insight  for  the  overall  control 
design,  but  it  has  also  provided  parameters  of  the  linearized 
dynamics . 

The  adaptive  controller  by  Schwartz  has  proved  to  yield 
satisfactory  behavior  in  terms  of  speed  of  response  and 
robustness  in  the  presence  of  model  changes.  On  the  other 
hand,  the  computational  complexity  involved  in  the  adaptive 
controller  calls  for  investigation  of  simpler  techniques  which 
do  not  require  complex  manipulation  of  matrices  to  be 
performed  in  real  time. 

With  these  motivations  in  mind  (mainly  adaptability  and 
robustness  in  the  presence  of  uncertainties) ,  we  have 
investigated  an  alternative  control  technique  based  on 
Variable  Structure  Control  (VSC) .  It  will  be  shown  that  a 
simple  and  robust  nonlinear  controller  can  be  designed, 
combined  with  an  adaptive  loop  which  compensates  for  some  of 
the  model  uncertainties. 

The  exciting  part  of  the  VSC  approach  is  that  the  design 
can  be  based  on  a  nominal  model  (usually  available  from  test 
runs  and/or  physical  insight)  and  a  bound  on  the  uncertainty 
of  the  model.  This  gives  to  the  designer  the  possibility  of 
using  this  simple  nominal  model,  possibly  linear  while  still 
preserving  the  global  stability  of  che  controlled  system. 
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Situations  which  exhibit  coupled  dynamics  can  be  handled  by 
VSC  design  techniques,  and  simple  and  robust  controllers  can 
be  designed  by  including  the  effects  of  coupling  in  the  model 
uncertainty. 

A  further  aspect  of  this  research  is  the  investigation  of 
the  performance  of  a  robust  linear  controller  in  the  presence 
of  uncertainties  and  nonlinearities.  In  particular  we  address 
the  problem  of  control  when  the  full  state  is  not  available 
for  measurement.  This  is  the  case,  for  example,  when  in  a 
diving  maneuver  we  measure  depth  alone,  rather  then  depth, 
pitch  and  pitch  rate.  It  turns  out  that  the  design  of  a 
"robust"  observer  (in  a  sense  which  will  be  clarified  in  a 
later  section)  together  with  an  optimal  controller  can  still 
be  satisfactory,  under  limited  ranges  of  the  operating 
conditions. 

An  outline  of  the  thesis  is  as  follows.  Chapter  II  covers 
modeling  of  the  AUV  and  stability  analysis,  while  Chapter  III 
addresses  the  controller  design  by  Variable  Structure  Control. 
In  Chapter  IV  a  controller  design  based  on  the  theory  of 
Robust  Observers  is  presented  with  applications  to  the  AUV. 
Alternative  hardware  realizations  are  surveyed  in  Chapter  V, 
while  conclusions  and  recommendations  for  future  research  are 
the  subjects  of  the  concluding  chapter. 
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II.  MODELING  OF  THE  AUV 


A.  BASIC  DEFINITIONS 

The  dynamics  of  underwater  vehicles  are  nonlinear  and  are 
affected  by  uncertainties  due  to  modeling  errors  (data 
sampling,  sensors,  electrical  noise,  physical  measurements, 
hydrodynamic  coefficients,  etc.)  and  external  disturbances 
(currents,  tether,  surfaces,  etc.).  In  general,  a  linear 
controller  for  a  vehicle  which  moves  in  "n"  degrees  of  freedom 
should  be  linearized  about  its  respective  axes  at  sufficiently 
different  speeds  and  then  some  form  of  "gain  scheduling"  is 

usually  implemented.  A  mathematical  model  of  the  AUV  can  be 

obtained  by  accounting  for  the  hydrodynamic  forces  acting  upon 
it.  In  particular,  the  model  results  in  a  set  of  highly 
nonlinear  differential  equations  in  state  space  form 

x  =  f(x,u,t)  (2.1a) 

where  x  is  the  vector  of  position  and  velocities  (both 
absolute  and  angular)  of  the  vehicle.  The  state  vector  x  is 
comprised  of 

x  =  surge  <f>  =  roll 

y  =  sway  6  =  pitch 

z  s  heave  ^  s  yaw 

which  represent  position  (absolute  and  angular)  in  a  global 
coordinate  frame,  and 
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u  =  surge  rate  p  s  roll  rate 

v  =  sway  rate  q  =  pitch  rate 

w  =  heave  rate  r  =  yaw  rate 

which  represent  velocities  in  a  body  fixed  coordinate  frame. 
Therefore  the  full  state  vector  x  in  the  differential  equation 
becomes  x=[u,v,w,p,q,r,x,y,z,^,0  ,t/>] 1 . 

In  a  diving  maneuver  the  only  states  of  interest  for  this 
diving  controller  are  depth  (z) ,  pitch  (©)  and  pitch  rate  (q) 
while  the  input  we  consider  is  the  dive  fin  (stern  plane) 
deflection  angle  (d) .  Therefore  the  dynamic  equations  become 

x=f (x, d)  (2.1b) 

where 

x=[q, ^ , z]  (2 . lc) 

and  d  is  the  stern  plane  dive  command.  Simulations  conducted 
using  the  full  nonlinear  model  [Ref.  2]  show  that,  within  the 
neighborhood  of  operating  conditions  (i.e,,  cruising  at  a 
constant  speed) ,  a  linear  model  yields  a  fairly  reasonable 
fit.  This  has  also  been  confirmed  by  data  collected  from  the 
prototype  [Ref.  1], 

The  basic  structure  of  the  dynamic  model  considered  is 
shown  in  Figure  2.  The  nonlinear  model  accounts  for 
perturbations  around  the  linear  model  which  are  caused  by 
various  hydrodynamic  effects  (currents,  cavitation,  surface 
effects,  etc.).  These  effects  have  the  tendency  to  perturb 
the  AUV  from  the  nominal  trajectory,  thus  acting  like  an  extra 
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Figure  2  AUV  Dynamic  Model 

component  on  the  stern  plane.  Therefore,  we  can  model  this 
effect  by  the  feedback  signal  fb2(t)  (shown  in  Figure  2)  which 
adds  a  perturbation  to  the  dive  plane  as 

fb2(t)  =  f(q(t),©(t),z(t),t).  (2.2) 
We  can  assume  that  the  perturbation  is  time-varying  and 
depends  on  pitch  rate,  pitch,  and  heave  (depth) . 

B.  LINEAR  MODEL 

To  develop  a  specific  linear  model  which  approximates  the 
dynamics,  let  us  consider  a  structure  of  the  form  shown  in 
Figure  3(a).  For  simplicity,  let  G(s)  be  a  first  order 
transfer  function  and  also  note  that  this  model  does  not 
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contain  a  feedback  from  0(t)  to  the  input  (i.e.,  the  gain  K0 
in  Figure  2) .  There  are  several  reasons  to  believe  that  this 
feedback  gain  is  required.  First,  physically,  the  equalizing 
moments  are 

Jq  =  F1(d(t))  +  F2(©(t))  +  F3  (2.3) 
where  F^dft))  is  the  torque  around  the  center  of  gravity  due 
to  the  stern  plane  angle,  F2(0(t))  is  the  torque  due  to  the 
displacement  between  the  Center  of  Gravity  (CG)  and  the  Center 
of  Buoyancy  (CB) ,  and  F3  is  the  balance  of  the  hydrodynamic 
moments  acting  over  the  body  of  the  vehicle  (including 
hydrodynamic  added  inertia,  drag,  and  angle  of  attack  induced 
moments) .  For  small  values  of  0(t) ,  this  effect  is 
proportional  to  0(t)  itself,  as  diagrammed  in  Figure  3(b). 
Second,  it  has  been  shown  in  Reference  2  that  the  response  of 
the  dynamics  relating  d(t)  and  0(t)  from  an  imparted  impulse 
to  the  stern  plane  yields  a  pitch  0(t)  which  decays  to  zero. 
This  indicates  that  all  the  poles  of  the  transfer  function 
between  d(t)  and  0(t)  are  stable.  Finally,  as  speed  increases, 
we  expect  that  the  stabilizing  effect  due  to  the  feedback  gain 
K0  decreases  because  other  hydrodynamic  forces  become  more 
dominant.  The  result  is  that  the  dominant  poles  of  'che 
transfer  function  from  d(t)  to  0(t)  follow  a  root  locus  path 
as  shown  in  Figure  3(c).  [Ref.  2] 

Following  this  formulation,  the  linearized  state-space 
mathematical  model  for  the  AUV  in  a  dive  maneuver  is 
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Figure  3  AUV  Characteristics 
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X 

= 
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(2 . 4c) 

where  we  define  the  "normalized  depth"  as 

z  =  ( **l/v)  z  (2 . 4d) 

with  v  being  the  forward  velocity  of  the  AUV.  The  fact  that 
the  system  is  in  controllable  canonical  form  will  be 
instrumental  in  the  design  of  the  control  system.  The 
parameters  a,  c  and  K0  in  (2.4a)  depend  on  the  operating 
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conditions  (speed,  primarily)  of  the  vehicle,  and  can  be 
determined  on  the  basis  of  test  runs  of  the  AUV.  On  the  basis 
of  these  runs  we  can  determine  nominal  values  a0,  c0  and  K0  for 
the  parameters  and  write  the  dynamics  as 

x  =  A0x  +  b(d  +f)  (2.5) 
and  design  the  controller  based  on  this  model.  In  the 
following  chapters  control  system  design  techniques  based  on 
this  model  will  be  presented. 

C.  STABILITY  METHODOLOGY 

Most  methods  of  stability  criteria,  including  Routh  and 
Nyquist,  are  not  applicable  to  nonlinear  systems.  The  second 
method  of  Lyapunov  is  the  most  general  test  for  stability 
analysis.  According  to  this  approach  we  define  stability  (in 
the  sense  of  Lyapunov)  and  a  criterion  for  stability  as 
follows. 

Definition:  Given  a  nonlinear  system 

X  =  f  (X,u,  t)  (2 . 6) 
and  an  equilibrium  point  xe,  such  that  f(xe,u,t)  =  0,  we  say 
that  xe  is  a  stable  equilibrium  point  (in  the  sense  of 
Lyapunov)  provided  that  for  all  e  >  0  there  exists  a  St  >  0 
(where  6  depends  on  e)  such  that  all  trajectories  of  (2.6)  for 
which  ||x(0)-xe||  <  e  do  not  leave  a  ball  of  radius  6C  around  xe, 
i.e.,  ||x(t)-xe ||  <  S(  for  all  t  >  0. 
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Definition:  If  xe  is  stable  and 

lim  x(t)  =  xe  (2.7) 
for  any  initial  condition  x(0) ,  then  xe  is  an  asymptotically 
stable  equilibrium  point.  [Ref.  3] 

A  method  to  determine  whether  an  equilibrium  point  is 
stable  or  not  is  by  the  second  method  of  Lyapunov.  This  method 
is  quite  convenient  for  the  stability  analysis  of  nonlinear 
systems  since  it  does  not  require  explicit  solutions  of  the 
differential  equations.  To  briefly  summarize  the  method, 
consider  a  system  described  by  the  differential  equation 
x=f(x,t)  with  f(0,t)=0  for  all  time.  If  there  exists  a 
differentiable  scalar  function  V(x,t)  such  that  V(x,t)  is 
positive  definite  and  V(x,t)  is  negative  definite,  then  the 
equilibrium  point  is  asymptotically  stable.  A  precise 
statement  of  the  method  is  presented  in  Reference  3 . 

Given  this  definition  of  stability  and  the  AUV  model 
presented  earlier,  we  were  ready  to  design  a  controller.  A 
nonlinear  design  method  using  Variable  Structure  Control  (VSC) 
was  then  followed. 
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III.  VARIABLE  STRUCTURE  CONTROL 


A.  INTRODUCTION 

The  nonlinear  dynamics  of  the  AUV  can  be  described  by 

x  =  f (x,u,t)  (3.1) 
and  present  a  high  degree  of  uncertainty  due  to  several 
factors.  In  spite  of  the  highly  nonlinear  behavior  we  can 
still  approximate  the  dynamics  of  the  AUV  by  a  linear 
component  and  a  nonlinear  perturbation 

x  =  Ax  +  b(d  +  f(x)).  (3.2) 
In  the  particular  case  of  a  diving  maneuver,  the  state  of 
equation  (3.2)  is  given  by  x  =  [q,0,F]‘  and  the  dynamics  (3.2) 
assume  a  controllable  canonical  form  structure  (see  (2.4a)). 

The  aim  of  the  controller  in  a  diving  maneuver  (like  in 
many  similar  problems)  is  to  drive  the  state  vector  x  to  track 
a  desired  state  xd.  In  other  words,  given  a  desired  state  xd 
e  of  the  form 

Xd(t)  =  [xd°(t)  ,♦  •  •  ,xd<rv1>(t)  ]  '  (3.3) 
we  want  to  design  a  controller  such  that  the  error  signal 

e(t)  =  x(t)  -  xd(t)  (3.4) 
tends  to  zero  as  t  -»  oo.  Furthermore  we  want  all  signals  in 
the  loop  to  be  bounded  for  all  time. 

In  our  particular  case  the  desired  signal  xd  is 
represented  by  the  desired  depth,  pitch  and  pitch  rate.  Since 
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in  steady  state  (i.e.,  cruising  at  a  constant  depth)  pitch  and 
pitch  rate  are  zero,  we  define  x<j  as 

xd  =  [0  0  -Zd/v]  '  (3.5) 

with  zd  being  desired  depth. 

Due  to  the  presence  of  the  uncertainty  f(x)  in  (3.2),  a 
control  solution  which  guarantees  a  sufficient  stabilty  margin 
has  to  be  devised.  The  Variable  Structure  Control  (VSC)  (also 
called  the  Sliding  Mode  (SM))  technique  is  at  the  basis  of  the 
controller  we  propose  for  the  AUV. 

B.  VSC  THEORY 

The  idea  behind  using  the  VSC  method  on  a  dynamic  system 

x  =  f (x,u,t)  (3.6) 

with  x  e  and  u  e  SR1  is  to  drive  the  state  from  any  initial 
condition  x(0)  onto  a  sliding  surface 

S (x)  =  0  (3.7) 

of  the  state  space  in  a  finite  time  and  to  keep  the  state  on 

the  surface  for  all  subsequent  times.  Although  the  surface 
(3.7)  can  be  chosen  in  a  fairly  arbitrary  fashion,  it 
nevertheless  must  be  associated  with  stable  dynamics.  In 

particular,  the  surface  s(x(t))  =  0  for  all  t  >  t0  must  imply 

lim  x(t)  =  0  (3.8) 

In  most  applications  the  sliding  surface  is  a  linear  function 
of  the  state,  in  the  sense  that  (3.7)  is  chosen  as 

s(x)  =  c'x  (3.9) 

for  some  vector  c*  e  9f. 
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From  the  statements  presented  in  the  last  paragraph,  the 
terminology  of  Eliding  Mode  is  evident,  since  the  state,  once 
taken  from  the  initial  condition  x(0)  onto  the  surface  s(x) 
=  0,  ’’slides"  to  zero  on  this  surface  by  virtue  of  (3.8)  (see 
Figure  4) .  Apart  from  the  requirement  of  being  associated  with 
stable  dynamics,  a  surface  in  is  a  sliding  surface  provided 


that  we  can  determine  a  control  input  signal  such  that 

.=  {x(t)  )s(x(t) }  <  0  (3.10) 

whenever  s(2f(t)  )  0.  The  reason  behind  condition  (3.10)  is 

that  the  definition  of  the  Lyapunov  function 

V(x)  =  0.5sz(x)  (3.11) 


yields 


along  the  trajectories  of  (3.6),  and  therefore  the  condition 

*  2 

V(x)  <  0  (given  by  (3.10))  makes  s  (x)  a  monotonically 

decreasing  function  so  that  s{x(t))  tends  to  zero.  Moreover, 
if  we  "strengthen”  condition  (3. 1C)  by  imposing 

s(x(t)}s{x(t)  )  <  -rj  |  s‘x(r.)  }  |  (3.13) 

with  r)  a  positive  constant,  then  it  is  easy  to  see  that  the 
sliding  surface  s(x)  =  0  is  reached  in  a  finite  time.  This 
can  be  seen  by  writing  (3.13)  as 

s{x(t)  }  <  -r?sgn{s (x(t)  }  (3.14) 

where  we  define  the  sgn  function  as 

+1  if  s  >  0 

sgn(s)  =  -1  if  s  <  0  (3.15) 

0  if  s  =  0 

On  the  basis  of  these  conditions,  if  we  restrict  ourselves 

to  linear  sliding  surfaces  (s(x)  =  c'x  as  in  (3.9)),  then  we 

obtain  from  (3.14)  and  (3.6) 

s{x(t)  }  =c'f(x,u,t)  <  -r??an{s(x)  }  (3.16) 

for  which  the  control  u(t)  must  be  sue.  ’at 

u+(t)  if  sgn{s(x  ,  >  0 
U (t)  =  (3.17) 

u.  (t)  if  sgn{s  (x)}  <  o 

and  u+  and  u.  are  such  that 

c'f(x,u+,t)  <  -r?  (3.18a) 

c'f(x,u.,t)  >  T]  (3.18b) 

with  r?  being  a  positive  constant  parameter  of  the  controller. 

As  a  matter  of  fact,  r?  does  not  have  to  be  a  constant  and,  in 

general ,  it  can  be  made  a  function  of  the  state  r/  (x)  .  From  the 
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definition  of  the  control  action  (3.17)  the  terminology  of 
Variable  Structure  Control  (VSC)  becomes  evident,  since  the 
control  assumes  different  structures  (u+  and  u.)  according  to 
which  side  of  the  sliding  surface  the  state  is  on  (i.e., 
according  to  the  sign  of  s(x)). 

C.  APPLICATION  OF  VSC  TO  TRACKING 

A  particularly  interesting  situation  arises  when  the 
dynamic  model  (3.6)  can  be  written  in  the  form 

x  =  f(x)  +  b(x)u(t)  (3.19) 
where  nominal  values  of  f  and  b  corresponding  to  f  and  b  in 
(3.,  19)  are  known.  For  simplicity  assume  that  b(x)  is 
completely  known  (this  will  be  relaxed  at  the  end  of  the 
section).  We  can  write  f(x)  as 

f(x)  =  f(x)  +  Af (x)  (3.20) 
with  Af  representing  the  uncertainty  in  f (x) ,  and  we  assume 
to  know  some  bound  on  | Af (x) |  to  be  specified  later.  Given  a 
desired  state  xd(t)  ,  we  want  to  apply  the  sliding  mode 
technique  in  order  to  move  the  state  x(t)  of  the  system  to 
track  the  desired  trajectory  xd(t)  .  In  order  to  do  this  we 
define  the  error  as 

e(t)  =  x(t)  -  xd(t)  (3.21) 
and  a  generic  surface  as 

s (e)  =  c'e  (3.22) 
for  some  vector  c'  e  sf.  Combining  (3.19),  (3.20),  (3.21)  and 
(3.22)  we  obtain 
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C1  (X  -  Xd)  =  c '  f  (x)  +  c'Af(x) 

-c'xd(t)  +c'b(x)u(t).  (3.23) 

An  important  feature  of  (3.23)  is  that  the  right  hand  side  is 
the  sum  of  terms  which  are  known  at  all  times  (c'f(x)  - 

c'Xd(t)  } .  an  uncertain  term  (c*  Af  (x) )  and  the  control  term 
(c'b(x)u(t) } .  For  this  reason  we  can  separate  the  control 
input  u(t)  into  two  terms  as 

U ( t )  =  U(t)  +  U(t)  (3.24) 

with  u(t)  compensating  for  the  "nominal"  dynamics,  i.e., 

c'f(x)  -  c'xd(t)  +  c'b(x)u(t)  =  0  (3.25) 

and  u(t)  compensating  for  the  uncertainty  Af.  Combining 
(3.23),  {->,.24)  and  (3.25),  we  obtain 

c ' e  =  c ' Af  +  c'bu  (3.26) 

and  we  can  determine  the  input  component  u  in  order  to  drive 
the  state  e(t)  on  to  the  sliding  surface  as 

s(e(t) }  =  c'e(t)  =  0.  (3.27) 

This  determination  can  be  accomplished  by  choosing  u(t)  so  as 
to  satisfy  (3.18),  which  leads  to  the  control 

u(t)  =  -F (x, t) sgn { c ' e (t) }  (3.28) 

with  F (x, t)  a  known  function  of  the  state,  such  that 

F (x, t)  >  |c'Af  (X)  |/|c»b(x)  | .  (3.29) 

Clearly  a  condition  which  has  to  be  satisfied  is  that 
c'b(x)  f  0  for  all  time. 
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In  conclusion,  if  c'  is  such  that  c*e(t)  =  0  (which 
implies  e(t)  ->  0) ,  then  we  can  say  that  the  control  input  from 
(3.28)  ,  (3.29)  and  (3.25) 

u(t)  =  u(t)  +  u(t)  =  - (1/c *b) (c1 f (x)  -  c'xd(t)  } 

+  F(x, t) sgn(c'e(t) )  (3.30) 

yields  the  following: 

•  the  error  e(t)  is  driven  to  the  sliding  surface 
c'e  =  C  in  a  finite  time,  and 

•  e(t)  -*  0  as  t  -•  w 

A  very  important  point  is  that,  once  the  state  is  on  the 
sliding  surface,  the  decay  of  e(t)  is  determined  by  the  vector 
c'  only.  [Ref.  4] 

D.  VSC  EXAMPLE 

To  illustrate  the  sliding  mode  design  process,  let  us 
implement  it  on  a  second  order  system  described  by 

x<2)  =  f(x,t)  +  b(x)u(t)  +  p(t)  (3.31) 

where  we  define  the  state  x(t)  =  [x(1)  x]  ' ,  the  desired  state 
Xd  =  [xd(1)  xd]  '  and  the  vector  c'  =  [1  A]  (with  A  >  0)  .  While 
f(x,t)  is  not  known  precisely,  we  assume  that  nominal  values 
f(x,t)  are  known,  for  which 

f (x,t)  =  f (x,t)  +  Af(x,t)  (3.32) 

with  Af(x,t)  having  a  known  upper  bound  of 

F  (X,  t)  >  |  Af  (x,  t )  |  .  (3.33) 

The  estimate  f(x,t)  is  available  from  several  sources  where, 
for  example,  a  nominal  model  could  be  generated  from 
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experimental  data  and  parameter  estimation  techniques.  If  the 
control  gain  b(x,t)  is  also  uncertain,  then  we  assume  that  it 
is  bounded  by  a  function  B(x,t)  where  the  gain  is  known  within 
a  certain  ratio 

1/B (x,  t)  <  b(x,t)/b(x,t)  <  B(X,t) .  (3.34) 

Also,  the  perturbations  p(t)  are  assumed  to  be  unknown  and 
bounded  by  a  continuous  time  function 

P(t)  >  |p(t) | .  (3.35) 

The  dynamics  of  s(x,t)  are  required  so  (3.9)  can  be 
differentiated  with  respect  to  time  to  get 

s(x,t)  =  x  -  xd  +  Ae.  (3.36a) 

Then  by  substituting  (3.31)  into  (3.36a)  we  obtain 

s(x,t)  =  f(x,t)  +  b(x,t)u(,t)  +  p(t)  -  xd  +  Ae  (3.36b) 
where  f(x,t)  and  b(x,t)  are  as  described  in  (3.32)  and  (3.34) 
respectively.  The  system  is  stable  if  s(x,t)  converges  to 
zero,  which  can  be  assured  if  u(x,t)  is  chosen  to  satisfy 
(3.13).  The  total  control  consists  of  two  parts  (as  in 
(3.24)),  one  for  the  known  or  estimated  part  of  the  dynamics 
(u)  and  the  other  for  the  uncertain  or  nonlinear  part  (u)  . 
Assuming  initially  that  there  is  no  uncertainty  in  either  the 
control  matrix  (b(x,t)  =  b(x,t))  or  the  transition  matrix 
(f(x,t)  =  f(x,t))  and  that  there  is  no  external  perturbation 
present  (p(t)  =  0),  then,  for  this  example,  the  known  control 
is  obtained  by  setting  (3.36b)  to  zero  so  that 

u(x,t)  =  - ( 1/b (x , t) )  [ f  (x, t)  +  Ae  -  xd] .  (3.37) 
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This  control  law  does  not  yet  satisfy  the  stability  criterion 
because  f  (x,t)5fcf  (x,t)  and  the  perturbations  p(t)  have  not  been 
accounted  for.  The  complete  control  law  will  have  another  term 
(u)  that  is  discontinuous  across  the  surface  in  order  to 
satisfy  ss<-»?|s|.  Now  with  all  the  uncertainties  included, 
where  s=s(X,t),  s  becomes 

s=f+Af+bu+p-xd+Ae  (3.38) 

and,  after  substituting  u=u+u  into  ss,  it  becomes 

ss  =  [Af  +  p  -  u]s  <  —77 1  s |  .  (3.39) 

This  results  in  u=-k(x, t) sgn(s)  so  that 

u(x,t)  =  u(x,t)  -  k(x,t)sgn(s)  (3.40) 

where  k(x,t)  is  determined  from  the  bounds  on  the 
uncertainties  and  perturbations  previously  estimated  as 

k(x,t)  =  [F(x, t)  +  P(t)  +  (3.41) 

This  variable  k(x,t)  is  the  mechanism  by  which  the  system 
uncertainties  are  accounted  for  and  which  will  cause  the 
discontinuous  part  of  the  control  to  compensate  for  their 
effect.  This  also  insures  that  s2(x,t)  is  a  Lyapunov  function 
which  in  turn  guarantees  stability. 

If  the  control  gain  is  also  uncertain,  then  the  control 
must  be  changed  to 

u(x,t)  =  [u(x,t)  -  k(x,t)sgn(s) ]/b(x,t)  (3.42) 

where  the  discontinuous  term  now  is 

k(x,t)  =  B [ F (x , t)  +  P(t)]  +  (B-l) |u(x,t) |  (3.43) 
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and  B  is  as  defined  in  (3.34)  [Ref.  5].  Again  notice  that  the 
control  discontinuity  has  become  larger  to  overcome  the 
increased  uncertainty  in  the  system. 

This  discontinuity  causes  a  chattering  problem  which  is 

present  due  to  the  abrupt  nature  of  the  sign  function,  so  a 

smoothing  operation  needs  to  be  implemented.  This  can  be  done 

by  replacing  the  sign  function  with  a  saturating  function 

which  has  a  user  defined  slope  or  boundary  layer  between  the 
« 

upper  and  lower  limits  of  plus  and  minus  one.  This  provides 
a  linear  region  around  the  surface  boundary  and  smooths  the 
control  action  (see  Figure  5a  and  5b) .  Now  the  designer  has 
another  variable  to  use,  <f>  (called  the  boundary  layex- 
thickness  and  whose  magnitude  depends  on  the  level  of  system 
uncertainties) ,  along  with  X  to  incorporate  into  his  design. 
These  variables  can  be  made  time-varying  to  account  for  times 
when  the  uncertainties  increase  or  decrease,  thus  increasing 
overall  performance.  [Ref.  6] 

The  questions  of  rigorous  mathematical  uniqueness, 
existence,  and  continuability  have  been  addressed  in  Reference 
5  and  Reference  7.  Though  bang-bang  controllers  are 
attractively  simple  and  optimal  control  specialists  have 
pretty  much  perfected  their  use,  the  problem  is  that  the 
differential  equations  governing  their  use  have 
discontinuities  on  the  right  hand  sides.  Therefore,  normal 
uniqueness  and  existence  theory  for  differential  equations  no 
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longer  applies.  Filippov  proves  the  existence  and 
continuability  for  his  own  solution  concept  [Ref.  7).  For  the 
question  of  uniqueness,  it  can  be  shown  that,  by  satisfying 
(3.13),  only  one  solution  to  the  discontinuous  differential 
equation  can  apply  at  any  given  time  [Ref.  5].  This  is 
equivalent  to  saying  that  the  derivative  of  the  state 
trajectory  must  always  point  towards  the  sliding  surface. 
Although  VSC  in  itself  is  a  powerful  design  methodology, 
incorporating  an  adaptive  portion  can  increase  the  robustness 
and  performance  of  the  controller. 
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E.  VSC  WITH  ADAPTIVE  COMPENSATION 


With  the  linear  model  developed  in  Chapter  II,  we  will 
now  determine  a  simple  adaptive  controller  capable  of  tracking 
the  state  of  a  reference  model  in  the  presence  of 
uncertainties  in  the  AUV  dynamic  parameters.  The  state  space 
description  of  this  model  is 

x(t)  =  Ax ( t)  +  bd(t)  +  f(x)  (3.44) 
where  d(t)  is  the  control  input  and  f(x)  is  the  nonlinearity 
which  comprises  model  inaccuracies,  external  perturbations 
and  changing  parameters.  This  depth  controller  design  is  based 
on  the  model  dynamics  with  an  adaptive  loop  to  place  the 
eigenvalues  for  the  desired  linear  part  and  a  switching  input 
based  on  VSC  design  for  the  nonlinear  part  f(x).  The 
controller  should  drive  the  AUV  to  a  desired  depth  zd  and  keep 
it  there,  despite  system  uncertainties  and  external 
disturbances.  The  difference  between  the  actual  state 
[q, 0, z/ (-v) ] 1  and  the  desired  state  [0, 0, Zj/  (-v)  ]  '  (where  v 
is  the  vehicle  speed)  is  defined  as  the  error  state 

e(t)  =  (q,e,  (z-zd)/(-v)  ] ' .  (3.45) 
Due  to  the  presence  of  integral  action  in  the  vehicle  model 
(3.44),  we  can  include  the  constant  desired  depth  in  the 
initial  condition  of  the  integrators  and  write  (3.44)  with  e 
as  the  state  (rather  than  x) . 

Now  let  us  choose  a  model  matrix  Aj„  with  eigenvalues 
determined  by  the  desired  closed  loop  response.  Clearly  Am 
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has  eigenvalues  in  the  stable  region.  If  we  choose  Am  such 
that  the  pair  (A^b)  is  in  controllable  canonical  form,  then 
it  is  a  simple  exercise  to  show  that  the  vehicle  dynamics 
(3.44)  can  be  written  as 

e(t)  =  A*e(t)  +  b (d (t)  +  K'e(t))  +  f(e)  (3.46) 
where  the  vehicle  dynamics  (at  the  current  operating 
conditions)  determine  the  required  gain  K.  This  dynamic  model 
(3.46)  is  at  the  basis  of  an  adaptive  controller  which  will 
drive  the  error  state  e(t)  to  zero  and  track  the  ordered 
depth.  The  surface  s(e)  =  0  is  defined  (as  described  in  (3.9)) 
as  a  linear  combination  of  the  error  states,  namely 

s(e)=c'e(t).  (3.47a) 
For  this  particular  formulation,  let  c'  be  the  left 
eigenvector  of  the  matrix  An,  associated  to  any  of  its  stable 
eigenvalues,  -A(i.e.,  c'An,  =  -Ac').  By  the  fact  that  Am  is  in 
companion  form  by  assumption,  it  is  possible  to  show  [Ref.  8] 
that  the  entries  of  the  vector  c'  are  the  coefficients  of  the 

polynomial  having  all  other  eigenvalues  of  Am  as  roots.  Since 

♦  • 

q=9,  9=z  and  zd  is  constant,  we  can  write  s(e)  as 

s(e)  =  czec2)  +  c^15  +  c0e  (3.47b) 
where  e=(z-zd)/ (-v) ,  and  c'  =  [c2  cy  c0] .  All  this  implies  that 
the  surface  s(e)=0  satisfies  one  of  the  requirements  of  being 
a  sliding  surface,  i.e., 

s  (e)  =0  **.’  lim  e  (t)  =  0  (3.47c) 
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On  the  basis  of  this  definition  we  can  determine  an  adaptive 
controller  which  drives  the  state  error  e  onto  the  surface 
s(e)=0  as  follows. 

In  the  model  (3.46)  let  us  assume  that  bounds 

Kjmin  <  Kj  <  KjMAX  (3.48) 

on  the  entries  of  the  vector  K,  and  a  bound 

F (e)  >  |c'f(e)/c'b|  (3.49) 

on  the  nonlinearities  are  known  to  the  designer.  Then,  under 
this  assumption,  the  control  input  is 

d(t)  -K'(t)e(t)  +  F(e)  sgn{s  (t)  }  (3.50) 

with  K  the  adaptive  gains  defined  as 

K(t)  =  -a (K(t)  }  -  /ie  (t)  s  (t)  .  (3.51) 

The  expression  aj{-K(t)}  is  given  by 

0  if  Kjmin  <  Kj  <  KjHAX 

aj(K(t))  =  -a  (Kj  (t )  -  Kjmin)  if  Kj  (t)  <  Kjrain 

-a  (Kj  (t)  -  KjHAX)  if  Kj  (t)  >  KjMAX  (3.52) 

(where  a  is  a  positive  constant),  hence  (3.51)  yields  a  closed 

loop  response  which  is  exponentially  stable  and  e(t)  -*■  0. 

[Ref.  4] 

Proof;  Using  the  fact  that  s=c'e  and  c'Am=-Ac',  we  can 
write  (3.46)  from  (3.47)  and  (3.50)  as 

s(t)  +  As (t)  =  c'bK' (t)e(t)  +  c'f(e) 

-c'bF(e) sgn{s(t) } ,  \  (3.53) 

where  K  =  K  -  K  is  the  parameter  error.  Define  the  Lyapunov 
function  as 

V(s,K)  =  0 . 5  (s2  +  glC'K)  (3.54) 
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where  g  =  c'b/p  (a  positive  quantity)  .  The  time  derivative  of 
(3.54)  is 

V(s,K)  =  -As2(t)  -  c'bF(e)  |  s(t)  |  +  c’f(e)s(t)) 

-gK'  (t)a(k(t)  ) .  (3.55) 

Now  we  can  see  that  the  term  K'  (t) a_ (K(t)  )  is  always 
nonnegative  by  definition  of  the  function  a(K(t)  = 
Oi  (K(t) )  ,• • •  ,ffn(K(t) ) ] '  in  (3.52).  Also,  by  the  definition 
of  the  bound  F(e)  on  the  nonlinearity,  the  term  (c'bF(e)  |s(t)  | 
-  c'f(e)s(t)}  is  always  nonnegative.  Therefore 

'  V(s,K)  <  -As2 (e  (t)  )  <  0  (3.56) 

along  the  trajectories  of  the  system.  This  implies  that 
s(e(t))  is  always  bounded  and  also  that  the  adaptive  gains 
K (t)  are  bounded.  This,  combined  with  (3.56),  yields 

Nf+Vdt  <  '  f+“  V (t) dt  <  V(0)-V(«»)  <  »  (3.57) 

0  •’o 

from  which  we  deduce  that 

lim  s(e(t))  =  0  (3.58) 

t-*o> 

and  e(t)  -»  0,  which  proves  the  result. 

F.  SIMULATIONS  AND  RESULTS 

The  controller  (3.50)  has  been  implemented  in  MATLAB  and 
the  basic  flow  chart  is  outlined  in  Figure  6.  The  performance 
is  satisfactory  at  all  speeds  and,  depending  on  the  desired 
closed  loop  eigenvalues  (i.a.,  the  choice  of  the  matrix  Am)  , 
the  response  could  be  tailored  towards  minimum  rise  time  or 
minimal  overshoot.  While  we  used  a  saturation  function  instead 
of  a  sign  function  in  our  program,  our  process  achieves 


26 


VSC  Algorithm 


Figure  6  VSC  Algorithm  Flow  Chart 

further  smoothing  of  the  discontinuous  boundary  by  use  of  the 
integral  process  in  determining  the  feedback  gains  Kj .  The 
simulation  results  are  shown  in  Figure  7  while  Figure  8  shows 
a  plot  of  typical  feedback  gains. 

These  results  show  that  the  controller  performs 
satisfactorily  over  a  wide  range  of  operating  conditions, 
ranging-  from  100  to  500  rpm  of  the  thrusters.  The  different 
rise  times  appearing  in  the  depth  plot  in  Figure  7  are 
consistent  with  the  different  speeds  of  the  vehicle.  In  all 
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Figure  7  VSC  Performance  Plots 

three  runs  the  dive  plane  action  (bounded  within  ±  0.4 
radians,  or  ±  23  degrees  due  to  the  physical  constraints)  does 
not  show  any  chattering  due  to  the  linear  saturation  adopted 
for  the  controller  (3.50),  rather  than  the  sgn  function 
formulated  in  the  theory.  The  result  is  a  smooth  operation  of 
the  dive  fin  (stern  plane) . 

The  adaptive  gains  Ki ,  K2  and  K3  plotted  in  Figure  8  stay 
within  reasonable  values  and  converge  to  constants  which 
depend  on  the  operating  conditions.  The  surface  S  is  also 
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plotted  and  is  shown  to  be  driven  to  zero  and  kept  there.  For 
these  particular  plots  the  speed  of  300  rpm  is  assumed. 

Considering  the  simplicity  of  the  controller,  we  could 
conclude  that  this  control  scheme  is  satifactory  for  on-line 
implementation  on  the  AUV  under  development  at  NPS.  Although 
the  VSC  method  is  an  efficient  nonlinear  design  method,  other 
linear  design  schemes  are  available  which  can  also  satisfy  the 
robustness  requirement. 
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IV.  LINEAR  ROBUST  CONTROL  OF  THE  AUV 


A.  INTRODUCTION 

The  VSC  described  in  the  previous  section  has  been  shown 
to  exhibit  satisfactory  robustness  properties  in  the  presence 
of  nonlinearities  and  unmodeled  dynamics.  A  major  requirement, 
however,  of  VSC  is  that  the  states  must  be  available  for 
measurement.  For  the  AUV  on  a  dive  maneuver  this  is  not  much 
of  a  drawback  since  the  state  signals  (pitch  rate,  pitch  and, 
depth)  are  provided  by  gyros  and  the  depth  cell. 

There  are  situations,  however,  in  which  we  might  want  to 
be  able  to  control  the  vehicle  with  incomplete  state 
information.  This  is  the  case  when  it  is  important  to  provide 
for  reliability  in  the  sense  of  ensuring  the  success  of  the 
mission,  even  in  the  presence  of  failure  of  a  gyro  or  its 
circuitry.  In  these  cases  a  model  based  on  an  observer 
provides  for  the  missing  measurements  by  estimating  them  from 
the  other  available  signals. 

It  is  a  well  known  fact  that  robust  design  techniques 
based  on  full  state  feedback  are  bound  to  lose  their 
robustness  properities  when  the  state  is  replaced  by  an 
estimate  from  an  observer.  This  has  been  pointed  out  by  Doyle- 
Stein  [Ref.  9]  in  the  context  of  robust  linear  controllers. 
In  the  next  section  we  introduce  the  notion  of  robust 
observers  designed  to  preserve  the  robustness  properties  of 
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the  compensator.  The  performance  of  the  AUV  with  t'.  is  type  of 
observer  is  given  in  the  last  section. 

B.  ROBUST  OBSERVERS  AND  THE  DOYLE-STEIN  CONDITION 

The  design  of  a  compensator  usually  starts  with  a  system 
description  such  as 

x  =  Ax  +  Bu  (4.1a) 

y  =  C'x  +  Du  (4 .  lb) 

where  A  is  the  plant  transition  matrix,  B  is  the  control 
matrix,  C  is  the  observation  matrix,  and  D  is  the  feedthrough 
matrix  (usually  it  is  zero) .  Because  of  the  Separation 
Principle,  the  controller  and  observer  can  be  designed 
independently.  Designing  the  controller  is  usually  the  first 
step  and  starts  with  finding  the  gains  of  the  feedback  control 
law  (u=-Gx)  where  the  required  response  is  determined  from 
system  specifications.  The  gains  are  then  selected  by  a  pole- 
placement  formula  such  as  the  Bass-Gura  method  or  Linear 
Quadratic  Regulator  (LQR)  methods,  based  on  a  quadratic 
performance  index,  and  hence  solve  a  Riccati  equation 
associated  with  the  system  (4.1)  .  Real  life  constraints  on  the 
physical  system  limitations  (such  as  the  power  supply)  limit 
the  size  of  the  control  u,  and  this  limitation  can  easily  be 
imbedded  in  an  LQR  design. 

When  all  state  signals  are  not  available  for  measurement, 
an  observer  provides  for  their  estimates  from  input/output 
measurements.  Using  a  Luenberger  observer  of  the  form 
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X  =  Ax  +  Bu  +  K(y-C'x)  (4.2) 
it  is  well  known  that  x~x  -»  0,  provided  A-Kc'  has  stable 
eigenvalues.  With  the  observer,  the  control  input  becomes 
u=-Gx.  The  use  of  an  estimated  state  (rather  then  the  state 
itself)  might  seriously  affect  the  robustness  of  the 
controller.  Even  the  seemingly  sensible  solution  of  making 
the  estimated  state  x  converge  to  the  actual  state  x  as  fast 
as  possible  (i.e.,  decrease  the  negative  real  part  of  the 
eigenvalues  of  A-Kc')  does  not,  in  general,  improve  robustness 
of  the  closed  loop  system  [Ref.  9]. 

A  significant  problem  concerning  a  lack  of  robustness 
(i.e.,  sufficiency  of  the  gain  and  phase  margins  in  the 
presence  of  parameter  variation)  has  been  associated  with 
systems  whose  observers  were  designed  solely  on  the  assumed 
sensor  noise  parameters  [Ref.  10].  Therefore,  the  issue  of 
robustness  of  the  closed  loop  dynamics  should  be  investigated 
as  part  of  the  design  process  of  an  observer.  A  typical  LQR 
controller  using  full  state  feedback  can  have  gain  and  phase 
margins  well  in  excess  of  six  decibels  and  60  degrees 
respectively;  however,  using  an  observer,  the  margins  can  be 
reduced  significantly  (refer  to  Figure  9)  .  The  effect  of  using 
faster  observer  poles  is  shown  in  Figure  9.  In  this  example 
(taken  from  [Ref.  9])  the  Nyquist  plot  of  the  system  with  full 
state  feedback  is  compared  with  the  plot  of  the  same 
controller  with  observed  state  feedback.  The  optimal  filter 
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(a  steady  state  Kalman  filter)  and  an  observer  with  "fast 
dynamics"  have  the  effect  of  greatly  diminishing  the  stability 
margins  of  the  system.  [Ref.  9] 


Figure  9  Nyquist  Plots  [From  Ref.  9] 

Consider  the  two  systems  in  Figure  10  where  both  systems 
are  minimum  phase,  observable  and  controllable.  For 
comparison,  the  controller  sections  have  been  highlighted  by 
placing  them  inside  the  dashed  lines.  Many  of  the  possible 
transfer  functions  have  been  compared  and  investigated,  with 
regard  to  robustness  from  the  connections  marked  X  and  XX.  It 
has  been  shown  that  the  closed  loop  transfer  functions  from 
r  to  x  are  identical  along  with  the  loop  transfer  functions 
(with  the  loops  broken  at  XX)  from  u*  to  u.  Yet  the  loop 
transfer  functions  (with  the  loops  broken  at  X)  from  u"  to  u* 
are  not  the  same  in  the  two  systems  unless  the  relationship 
K[I  +  C(SI  -  A)  ^K]  1  =  B[C  (SI  -  A)  ’'b]  *1  (4.3) 
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Figure  10  Feedback  (a) Full-State  (b) Observer  [After  Ref.  9] 
holds.  Recall  that  K  is  the  observer  gain  which  now  has 
certain  restrictions  placed  on  it  from  (4.3).  This 
relationship  is  known  as  the  Doyle-Stein  (D-S)  condition,  is 
independent  of  the  control  gain  and  depends  only  on  the  open- 
loop  characteristics  of  the  observer.  When  (4.3)  holds,  the 
observer  does  not  influence  the  transfer  function  from  r  to 
x  [Ref.  10],  Finally,  note  that  the  loop  transfer  functions 
are  the  same  only  at  point  XX,  which  is  inside  the  designed 
controller,  but  they  are  not  the  same  at  point  X,  which  is  the 
control  interface  to  the  plant  and  where  natural  parameter 
uncertainties  will  most  probably  occur.  [Ref.  9] 

Doyle  and  Stein  have  shown  that,  if  the  observer  gains  are 
parameterized  as  a  function  of  q  (a  scalar  variable) ,  then 


this  particular  function  K(q)  must  meet  certain  requirements 
as  q-oo  with  the  main  one  being  that 

K(q)/q  =>  BZ  (4.4) 
where  Z  is  any  non-singular  matrix  [Ref.  9).  A  Kalman  filter 
provides  for  such  a  function  (one  of  many  possibilities)  that 
works  for  all  controllable,  observable  and  minimum  phase 
systems.  This  is  similar  to  assuming  that  the  covariance  Q  of 
the  process  noise  used  in  the  Kalman  filter  has  the  special 
form  of 

Q(q)  =  Qo  +  q2BWB'  (4.5) 
where  W  is  any  positive  definite  symmetric  matrix  and  q  is 
the  scalar  weighting  factor  defined  above.  (The  value  q=0  is 
associated  with  the  nominal  steady  state  Kalman  filter  gain 
for  the  assumed  process  noise  Q0. )  As  q-oo,  the  response  tends 
toward  the  full  state  feedback  response  as  shown  in  Figure  11. 
The  extra  term  in  (4.5)  is  equivalent  to  adding  extra  process 
noise  to  the  control  input  of  the  plant.  Normally  q  is 
increased  until  a  satisfactory  compromise  between  robustness 
and  noise  performance  is  obtained.  This  method  gives  an  easy 
way  to  balance  the  requirement  between  noise  rejection  and 
stability  margins.  We  have  investigated  this  method  with 
regard  to  the  NPS  AUV  vehicle  model  and  present  the  results 
in  the  next  section.  [Ref.  9] 
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Figure  11  Fictitious  Noise  Nyquist  Plots  [From  Ref.  9] 

C.  SIMULATIONS  AND  RESULTS 

Following  the  background  development  of  a  Doyle-Stein 
observer  (robust  observer)  that  was  presented  in  Chapter  II, 
we  first  need  to  generate  numerical  values  for  the  linear 
model  of  the  AUV.  Presently,  there  are  many  simulation 
languages  available,  but  we  used  MATLAB.  The  complex  and  large 
C  program  that  encompasses  the  AUV  model  was  implemented  as 
a  function  in  MATLAB.  This  allowed  us  to  try  many  forms  and 
types  of  controllers  and  to  simulate  them  in  a  much  easier 
fashion  than  a  conventional  language  might  have  afforded. 
Mainly  this  was  because  of  the  rich  libraries  and  simple 
graphics  that  MATLAB  has  and  the  relative  ease  of  changing  and 
modifying  the  algorithms  of  the  controller. 

The  AUV  model  can  be  fit  to  an  A.RX  (AutoRegressive) 
discrete  time  model  of 
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Atq'^yCkT)  =  BCq'^ufkT)  +  e(kT)  (4.6)' 

where  A  and  B  are  polynomials  in  the  time  delay  operator  q*1, 
T  is  the  sampling  interval  and  e  is  an  error  sequence 
(possibly  colored)  .  Using  this  ARX  model  and  a  Recursive  Least 
Squares  (RLS)  or  Recursive  Instrumental  Variable  (RIV) 
algorithm,  the  parameters  of  the  model  can  be  estimated.  In 
general  the  RIV  methods  are  better  suited  to  cases  when  the 
noise  is  significantly  colored  (non-white) .  In  our  case,  a 
third-order  model  has  been  estimated  using  both  RLS  and  RIV 
with  only  small  differences  noted  in  the  numerical  results. 
The  input  to  the  RLS  algorithm  is  required  to  be  "persistently 
exciting"  so  that  all  modes  of  the  model  will  be  excited. 

Since  a  fourth  order  model  yields  very  similar  results, 
the  choice  went  to  the  simpler  third  order  model,  further 
confirming  the  model  choice  made  in  Chapter  II.  Figure  12  is 
a  plot  of  the  input,  actual  output  and  estimated  output. 

The  next  step  was  to  generate  the  state  feedback  and 
estimator  gains.  MATLAB  functions  were  used  extensively  in 
this  area.  For  state  feedback  any  pole  placement  techniques 
will  work  such  as  DLQR  or  PLACE  in  MATLAB.  We  used  the  steady 
state  Kalman  gains  from  DLQE  for  the  estimator  gains  where  the 
process  noise  was  adjusted  or  modified  according  to  equation 
(4.5).  Copies  of  the  MATLAB  programs  used  in  this  thesis  can 
be  obtained  by  contacting  Professor  Roberto  Cristi  (fourth 
entry  in  the  initial  distribution  list,  phone  408-646-2223) . 
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Figure  12  RLS  Model  Input/Output  Characteristics 
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Figure  13  shows  the  performance  of  the  controller  for 
several  choices  of  robust  observer.  As  outlined  in  the 
previous  section,  the  robust  design  yields  a  family  of 
observers  (parameterized  by  the  parameter  q)  which  improve 
robustness  as  q  becomes  larger.  This  effect  is  shown  in  Figure 
13,  where  the  depth  responses  obtained  with  q=0,  10  and  50  are 
shown.  Notice  the  increased  stability  and  reduction  of 
oscillations.  This  is  also  shown  in  the  input  signal  which  is 
limited  to  ±  0.4  radians. 

Regardless  of  the  type  of  controller  used  (VSC  or  robust 
observer) ,  an  efficient  implementation  in  hardware  requires 
at  least  a  survey  of  the  available  AUV  architectures. 
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Figure  13  D-S  Observer  Performance  Curves 


IV.  AUV  ARCHITECTURE 


A.  INTRODUCTION 

AUV  systems  and  their  related  hardware  have  become 
increasingly  complex  in  order  to  satisfy  all  the  levels  of 
vehicle  control.  These  levels  encompass  varying  degrees  of 
asbtraction  from  the  highest  (Artificial  Intelligence)  to  the 
lowest  (actuator  control  algorithms) .  Historically,  the 
challenge  has  been  to  find  some  sort  of  hardware  and  software 
combination  to  satisfy  all  the  constraints  generated  by  such 
a  sophisticated  system.  The  problem  is  that  a  general  purpose 
microcomputer  solution  tends  to  be  slow  (non-real  time)  and 
inefficient  or  non-optimal.  The  trend  is  to  "tailor"  the 
hardware  and  software  to  the  problem  at  hand.  This 
"application  specific"  approach  is  especially  germane  with 
robotics  and  AUV  systems  since  they  rely  heavily  on  sensory 
feedback  where  real-time  response  is  a  necessity.  Most  of  the 
systems  surveyed  were  set  up  in  some  sort  of  hierarchical 
fashion  corresponding  to  the  various  levels  of  "intelligence" 
or  abstractions  in  their  mission  plans.  Each  level  has  its 
own  numerical  computational  requirement;  therefore,  hardware 
and  software  selections  for  each  level  must  be  tailored  for 
that  specific  application.  This  should  allow  the  system  to 
be  as  efficient  as  possible  and  to  allow  the  lowest  control 
level  to  run  in  real-time.  This  chapter  looks  at  some  of  the 
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various  hardware  methods  which  can  help  achieve  this  real-time 
performance  requirement  at  the  control-actuator  interface  and 
yet  still  allow  some  flexibility  in  design  which  will 
encompass  a  wide  range  of  applications. 

A.  HARDWARE  ALTERNATIVES 

The  architecture  specifications  should  match  or  exceed 
the  performance  requirements  of  the  algorithm  which  are 
normally  formulated  in  terms  of  latency  (elapsed  time  from 
when  the  input  is  present  until  the  output  is  ready)  and 
throughput  rate  (given  in  MIPS  or  MFLOPS) .  There  is  a 
tendency  to  concentrate  on  the  throughput  performance  while 
the  latency  (i.e,  delay)  is  of  critical  importance  in  control 
algorithms  (the  system  stability  requirements  may  allow  only 
so  much  inherent  delay) .  Control  applications  normally  require 
positional  accuracy,  concurrent  tasks,  repeatability, 
robustness  and  timing  constraints  (tasks  must  all  be  done 
within  the  sampling  period)  for  the  hardware.  It  seems 
reasonable  that  the  optimal  approach  would  be  some  combination 
of  the  various  hardware  technologies.  Some  of  these 
technologies  are: 

•  pipelining, 

•  RISC  (Reduced-Instruction-Set-Computer) , 

•  bit-sliced  microprocessors, 

•  vector  processing, 

•  DSP  (Digitial  Signal  Processing)  chips, 
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•  ASIC  (Application-Specific  Integrated  Circuit)  chips, 

*  systolic  arrays, 

♦  multiple  processors,  and 

•  PLDs  (Programmable  Logic  Devices) .  [Ref.  11] 

Pipelining  (paralleling  the  datapath)  tries  to  improve 

the  throughput  by  shortening  the  clock  cycle,  but  then  the 
latency  also  increases.  The  concept  will  also  work  well  at 
higher  levels  of  control  or  abstraction  (concurrent 
programming)  and  not  just  at  the  computational  level  such  as 
a  multiplier/accumulator  (MAC) .  [Ref.  ll] 

RISC  is  a  diverse  technology  and  is  still  subject  to  much 
debate  and  company  specific  design  philosophy.  The 
performance  comes  from  using  less  "chip  real  estate"  to  encode 
fewer  instructions  thus  leaving  more  room  to  add  additional 
components  to  speed  up  overall  program  execution.  Several 
manufacturers  have  working  RISC  microprocessors  and  their 
throughput  performance  is  impressive.  RISC  systems  have  not 
had  many  commercial  applications  because  twice  as  much  memory 
(versus  a  80X86  system)  is  required  and  their  speed  advantage 
has  eroded  since  the  newer  80X86  chips  are  much  faster  now. 
Even  so,  AUV  control  is  an  area  where  a  dedicated  RISC 
architecture  might  be  exploited  in  a  generic  controller 
scenario. 

Bit-sliced  microprocessors  are  very  fast  since  they 
generally  use  an  ECL  (Emitter-Coupled-Logic)  chip  set.  They 
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can  be  tailored  to  an  individual  application  and  fit  very  well 
into  fault  tolerant  systems.  Their  power  hungry  ECL 
construction  and  multi-chip  expense  limit  their  general  use 
in  a  space  and  energy  conservative  environment  such  as  an  AUV, 
but  certain  performance  constraints  requiring  ultra  fast 
processing  would  benefit  from  their  use.  [Ref.  12] 

Systolic  arrays  and  vector  processing  basically  combine 
several  arithmetic  or  processing  units  in  various  geometries 
to  improve  the  data  stream  flow.  The  general  idea  is  to 
increase  performance  by  some  sort  of  parallel  computation. 
Most  of  the  supercomputers  use  some  form  of  vector  processing 
to  achieve  enormous  throughput  capability  and  systolic  arrays 
have  been  used  fairly  successfully  in  image  processing 
systems.  The  multi-dimensional  nature  of  an  AUV  system  along 
with  the  MIMO  control  problem  (with  its  discrete  time 
dynamical  difference  equations)  requires  a  matrix  formulation 
and  then  a  real-time  solution  which  lends  itself  to  this  type 
of  processing. 

DSP  chips  have  begun  to  flourish  and  their  performance  is 
truly  remarkable.  Several  of  the  well-established  companies 
that  market  32-bit  floating  point  processors  are:  Motorola 
(DSP96002),  Texas  Instruments  (TMS320C30)  and  AT&T  (DSP32C- 
80)  .  While  the  state-of-the-art  DSP  chips  are  relatively 
expensive,  their  performance  (33  MFLOPS  for  the  T.I.  chip), 
onboard  memory  and  on  chip  I/O  make  them  virtually  a  single 
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chip  solution  for  almost  any  control  or  filter  problem.  The 
capability  of  extensive  on-chip  programming  and  the  ability 
to  move  blocks  of  memory  onto  the  chip  greatly  enhances  the 
actual  throughput  of  a  large  control  algorithm.  A  possible 
pitfall  with  these  chips  is  that  along  with  the  substantial 
performance  increases  come  commensurate  leaps  in  system 
details  (debugging  programs,  etc.).  The  most  effective  way 
to  minimize  this  problem  and  its  learning  curve  is  to  purchase 
the  complete  package  (C  compiler,  PC-board,  etc.)  from  the 
manufacturer.  Here  again  it  is  probably  best  to  stay  with 
well-established  companies  who  are  marketing  their  next 
generation  chip  which  already  has  most  of  the  required 
software  available  from  the  previous  generation's  chip  product 
line.  The  final  choice  of  DSP  chip  will  be  a  balance  between 
performance,  cost  and  supporting  software  (which  is  the  most 
important  item  considering  the  essentially  equivalent 
performance  from  all  of  the  32-bit  chips)  .  The  time  and 
effort  involved  in  getting  any  DSP  chip  to  actually  work  in 
a  system  will  entail  a  large  portion  of  the  total  resources, 
so  any  product  that  can  shorten  or  automate  this  process  will 
be  money  well  spent. 

The  term  ASIC  seems  to  suggest  that  it  is  a  combination 
of  some  design  methodology,  system  requirements  and  choice  of 
algorithm  which,  taken  collectively,  purport  to  solve  some 
specific  control  problem.  This  would  allow  many  of  the  system 
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overhead  functions  usually  performed  by  a  general  purpose 
microprocessor  to  be  eliminated,  and  hence  an  overall 
performance  increase  would  be  achieved.  ASIC  devices 
typically  are  either  full  custom  or  semi-custom  with  the 
latter  the  most  common  due  to  the  greatly  reduced  design 
costs.  Semi-custom  designs  normally  use  gate  arrays  or 
standard  cells  (gates,  registers,  ALU,  etc.)  with  the  former 
the  most  popular  since  it  generally  has  a  faster  design  turn 
around  time.  A  natural  application  for  these  chips  would  be 
in  an  AUV  system  where  navigation,  actuator  controllers  and 
sensor  processing  could  require  vast  computational  power. 
Generally  the  problems  to  be  addressed  are  testability  of  the 
ASIC  design  and  efficient  layout  (interconnections)  of  the 
chip  to  minimize  propagation  delays.  [Ref.  13] 

One  of  the  more  exciting  areas  to  explore  is  that  of 
multiple  processors.  There  are  various  methods  and  ideas  on 
how  to  implement  these  distributed  computing  techniques.  One 
of  the  more  important  decisions  is  the  one  between  message 
passing  and  shared  memory.  Message  passing  can  have  a  long 
latency  time  and  a  large  overhead  due  to  protocol 
requirements.  The  allure  of  multiple  processors  is  that  the 
type  of  processor  used  can  be  chosen  to  best  complement  the 
particular  algorithm  at  hand.  How  the  inter-communications 
between  processors  is  handled  is  a  crucial  design  decision 
second  only  to  the  processor  choice  itself.  [Ref.  11] 
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There  are  essentially  only  a  few  methods  of  inter¬ 
communication  between  processors:  RS232  cable,  a  local  or 
system  bus,  shared  memory  (dual  port  RAM,  etc.)  or  a  local 
area  network  (LAN)  line  (e.g.,  E\  srnet) .  The  LAN  tends  to 
have  a  lot  of  overhead,  is  serial  in  nature  and,  hence,  is 
relatively  slow.  For  distant  processors  that  need  infrequent 
communications  this  could  be  a  viable  method,  but  for  tightly- 
coupled  processors  a  versatile  bus  is  the  only  real  choice. 
There  are  several  styles  and  types  of  buses  from  which  to 
choose;  they  range  from  simple  passive  backplanes  to  forward- 
looking  high-speed  32-bit  buses.  Issues  such  as  compatibility 
and  performance  must  be  investigated  and  the  inevitable 
compromises  made.  Several  types  of  buses  can  be  employed  in 
one  system  since  AUV  control  systems  tend  to  be  hierarchical 
in  topology  and  each  layer  of  control  or  abstraction  can  be 
tied  together  with  a  particular  bus  to  fulfill  a  specific 
need. 

The  newest  32-bit  buses  offer  the  most  performance  and 
versatility  so  far.  They  can  block  transfer  to  RAM,  do  cache 
coherence,  autoconfigure  (poll  boards  attached  to  the  bus  and 
then  adjust  the  related  interface  software)  and  interact  with 
the  fault-tolerant  system  (i.e.,  logically  remove  the  faulted 
board  on-line)  .  Some  of  the  buses  can  handle  larger  bus 
widths  of  up  to  256  bits  of  data  (e.g.,  FUTUREBUS+  :  IEEE 
standard  P896)  .  The  U.S.  Navy  has  decided  to  base  all 
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mission-critical  computers  on  FUTUREBUS.  Given  this  decision 
and  the  superb  performance  and  adaptability  of  this  bus,  it 
would  seem  a  natural  choice  for  future  expansions.  Of  special 
note  is  the  new  type  of  transceiver  used  by  FUTUREBUS  called 
BTL  (Backplane  Transceiver  Logic) .  This  transceiver  reduces 
the  bus  capacitive  load  and  hence  increases  the  bus  bandwidth 
to  400  Mbytes/sec  (an  order  of  magnitude  better  than  ECL 
transceivers) .  [Ref.  14] 

Finally,  PLDs  of  which  the  EPROMs  and  PALs  are  examples 
are  a  relatively  cheap  way  to  implement  functions  or 
processes.  PLDs  are  expandable,  universally  compatible  and, 
by  design,  tailor-made  to  the  specific  application  and 
algorithm.  These  devices  are  cheap  (compared  to  DSP  devices) 
and  several  versions  are  in  development  that  have  provisions 
to  be  re-programmed  on  the  fly  (EEPROMs,  etc.). 

B.  OPERATING  EXAMPLES 

While  there  are  as  many  theories  and  ideas  on  how  an  AUV 
architecture  should  be  built  as  there  are  AUV  manufactorers, 
it  can  be  quite  helpful  to  investigate  a  few  designs  which 
have  been  built  and  are  actually  operational.  This  section 
will  briefly  overview  three  vehicles: 

•  ARCS  (Autonomous  Remotely  Controlled  Submersible)  built 
by  International  Submarine  Engineering  Limited  (ISE)  in 
Canada. 

•  EAVE  (Experimental  Autonomous  Vehicle)  EAST  built  by  the 
University  of  New  Hampshire. 
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•  FS  (Free-Swimmer  and  formally  EAVE  WEST)  built  by  Naval 

Ocean  Systems  Center  in  San  Diego,  CA. 

Table  1  is  a  pc\rtial  list  of  several  AUV  designs  and  their 
ROV  manufacturers.  Two  comprehensive  AUV  references  that 
detail  currant  vehicles  are  ROV  Review  1990  and  Undersea 
Vehicles  Directly  1990.  (ROV  Review  can  be  purchased  after 
January  1990  by  contacting  Dean  Given,  PO  Box  368,  Spring 
Valley,  CA  92077,  phone  619-660-0402.  Undersea  Vehicles 
Directory  can  be  purchased  after  November  1989  by  contacting 
Frank  Busby  in  Arlington,  VA,  phone  704-892-2888.)  [Ref.  15] 

The  ARCS  is  a  commercially  available  AUV,  so  it  uses  "off- 
the-shelf"  hardware  and  software  .to  reduce  costs.  The 
modeling  of  the  control  system  organization  followed  that  of 
a  typical  naval  submarine  and  a  multi-tasking  setup  was  used 
to  schedule  all  the  tasks.  The  specific  hardware  included  a 
16-bit  CPU  for  the  multi-tasking,  an  Intel  Multibus  and  three 
single  board  computers  (two  8086s  with  8087s  and  one  8088), 
tied  together  with  a  common  backplane  bus  and  sharing  some 
dual  port  RAM.  Interfaces  to  external  equipment  were  through 
RS-232  serial  links.  The  requirements  for  a  real-time  multi¬ 
tasking  operating  system  led  ISE  to  choose  Intel's  RMX86 
system.  [Refs  16  and  17] 

The  University  of  New  Hampshire  has  been  developing  and 
enhancing  their  AUV  project  since  1977. __  They  .are  currently 
pursuing  their  third  generation  vehicles  which  stress 
Artificial  Intelligence  (A. I.)  concepts  and  multi-vehicle 
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Table  I  Current  AUVs  [From  Ref.  15] 


I  wn*  1 

cn*T1 

turn 

Wrtncad  Under**# 

*  \*h»cto  (A.UVJ 

KA 
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Gad  el  Mafco  dewonstnmg  high  mjnet/vetabCtty.  but  none  since  September  1963 
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CarneTe-Menon  University 
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62.5  mmee  project  underway  lo  devttop  vehide  by  1967  *nth  real  Dme  xoustlc  speed  erf 

2  knots,  ter  Inspection  use  on  surveys  and  around  offshore  oil  pfatlorms 

Erperknenui  Autonomous 
VlMcU  £jii  (lave  Cist) 

900 

University  0 (  New  Hampshire, 

Our  ham,  N  H 

A  1.5-mekf  bo*  tee  autonomous  vtWde  test-bed  lot  structure  and  otpttlni  Inspection: 
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mous  cowrol  with  te*emeiry  override:  unmanned  veMdu  mainly  lor  prootrfsion  stupes  but 
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Hwloi'Wjii  University. 
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University  d  Washington. 
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U  S  Naval  Research  laboratory. 
Washington.  0  C 

Pidotype  if  program  wttated  in  1975  to  dnetop  Ngh  speed  seafrh  vthtcie  with  rang* 
or  1900»*.  has  undergone  shadow  water  tests;  project  currenity  inactive 
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1000 
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Marseilles,  franca 

Prototype  ter  structural  Inspection  and  drlling  support:  several  tests,  program  currently 
Inaclhe 
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teoiolype  designed  for  noduk  collection 

cooperation.  EAVE  EAST  III  has  much  greater  memory  capability 
and_ processing  power  than  its  predecessors  in  order  to  serve 
as  a  testbed  for  future  concepts  and  evaluations.  Their 
computer  system  is  hierarchical  with  three  distinct  levels. 
The  lowest  level  handles  all  the  fast  (less  than  one  second) 
systems  including  reading  sensors,  pre-processing  data  and 
controlling  actuators.  This  level  consists  of  three  68000 
processors  which  communicate  over  RS-232  cable  to  each  other 
and  to  an  interface  with  the  higher  levels.  This  lower  level 
can  control  the  vehicle  without  any  higher  level  assistance. 

The  higher  level  architecture  consists  of  several  68020 
processors  on  a  VME  bus  and  liberal  *use  of  dual  port  memories 
which  allow  local  CPU  use  as  well  as  communication  by  other 
CPUs  on  the  bus.  Figure  14  shows  this  multi-bus  setup  which 
helps  to  incorporate  I/O  and  data  storage  into  the  overall 
system.  Their  choice  of  a  real-time  operating  system  was  pSOS 
which  was  also  tasked  with  running  a  symbolic  language  (LISP)  . 
So  far  their  choice  of  architecture,  hierarchy  and  software 
has  proven  to  be  reliable,  capable  and  extensible.  [Ref.  18] 

NOSC  has  a  wealth  of  experience  in  underwater  vehicles 
and  their  subsystems.  Their  FS  (Free-Swimmer  or  previously 
known  as  EAVE  WEST)  vehicle  has  been  through  several  upgrades 
in  capability.  Although  it  is  probably  not  as  advanced  as 
some  of  their  current  vehicles  (which  have  restricted 
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Figure  14  EAVE  EAST  Architecture  [From  Ref.  18] 

fairly  representative  of  their  architectural  design 
philosophy.  The  NOSC  undersea  branch  is  committed  to  all 
aspects  of  AUV-related  research  from  using  a  transputer  array 
to  simulate  an  Artificial  Neural  Network  (ANN)  to  designing 
and  testing  a  Plan  Execution  System  (PES) .  [Ref.  19] 

The  FS  architecture  takes  a  layered  approach  as  pictured 
in  Figure  15.  Their  PES  runs  on  an  80286  CPU  and  communicates 
with  the  navigator  CPU  (8088)  and  the  lower  level  bus  through 
a  RS-232  link.  The  lower  level  groups  several  functions 
(using  8088  CPUs)  together  along  with  sensors  and  actuator 
(effector)  I/O.  Their  system  requirements  lead  them  to  use 
VERTX  for  the  operating  system. 
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Figure  15  NOSC  FS  Architecture 


This  is  just  a  sampling  of  the  current  operational  AUV 
architectural  designs.  Table  2  lists  several  other 
organizations  and  some  of  their  research.  [Ref.  12]. 


D.  OBSERVATIONS 

The  selection  of  an  AUV  architecture  is  a  very  critical 
process  because  the  system  must  be  expandable  and  compatible 
enough  to  meet  ever-increasing  complex  missions  needs.  The 
design  of  a  generic  AUV  architecture  would  entail  several  key 
items . 

First,  the  overall  system  would  be  hierarchical  in  nature 
with  several  layers  of  control.  This  aspect  would  suggest 
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Table  II  AUV  Research  Organizations  [After  Ref.  12] 


VISITED 


PlIRPflSiF/ IMPORTANCE 


AMETEK-STRAZA/SAN  DIEGO 
ARL/UNIV.  OF  TEXAS/AUSTIN 
AUT0NET1CS,  LOS  ANGELES 
DRAPER  LABORATORIES,  BOSTON 


HYDROPRODUCTS,  SAN  DIEGO 

EDO  WESTERN.  SALT  LAKE  CITY 
GOULD.  BALTIMORE 

HONEYWELL  MARINE,  SEATTLE 

LOCKHEED  ADVANCED  MARINE 
SYSTEMS,  SUNNYVALE 

MARTIN  MARIETTA,  BALTIMORE 
MARTIN  MARIETTA,  DENVER 
NOSC,  SAN  DIEGO 
NOSC.  HAWAII 


PERRY  OFFSHORE.  RIVIERA  BEACH 

UNIVERSITY  OF  UTAH;  SALT  LAKE  CITY 
WESTIHGIIOUSE,  ANNAPOLIS 

DEEP  SUBMERGENCE  LABORATORY  OF 
WOODS  HOLE  OCEANOGRAPHIC 
INSTITUTE.  FALMOUTH 


PROVEN  CAPABILITIES,  TRENDS.  WORK  SUITES 
OBSTACLE  AVOIDANCE  SONAR  ; 

NAVIGATION,  AOVs,  HYDRODYNAMICS,  LOW  DRAG  BODIES 

AUTONOMOUS  VEHICLE  PROCESSING  AND  CONTROL, 
NAVIGATION,  RELIABILITY  ANALYSIS,  LIGHTWEIGHT 
STRUCTURES 

PROVEN  CAPABILITIES.  TRENDS,  UNTETHERED 
APPLICATIONS 

SONARS,  TV  CAMERAS 

ENERGY  SUPPLIES.  FIBER  OPTICS,  CABLE  HANDLING 
SYSTEMS 

NAVIGATION,  ACOUSTIC  LINKS,  VEHICLE  CAPABILITIES, 
TRENCHING,  MK-50  TORPEDO.  BURIED  OBJECT  DETECTION 

LOW  COST  MINE  NEUTRALIZATION  VEHICLE 


LIGHTWEIGHT  STRUCTURES,  BONDING,  AOVs,  lRRD 
PROGRAMS 

INTELLIGENCE  PROCESSING,  ROBOTICS.  AUTONOMOUS  LAND 
VEHICLE 

UNMANNED  VEHICLE  OVERVIEW.  AOVs,  ACOUSTIC  LINKS, 
WORK  FUNCTIONS 

DEEP  VEHICLES.  AUTONOMOUS  VEHICLES,  FIBER  OPTIC 
CABLES.  CABLE  DEPLOYMENT  METHODS,  MARINE 
BIOSYSTEMS 

COMMERCIAL  VEHICLES,  MISSION  SPECIFIC  VEHICLES, 
WORK  PACKAGES,  MANIPULATORS 

DEXTEROUS  HAND 

PROVEN  CAPABILITIES,  LAUNCH/RECOVERY  SYSTEMS, 
SUBMERSIBLE  VEHICLE  INTERFACE,  SONARS 

DEEP  SUBMERGENCE  VEHICLES,  IMAGING  SYSTEMS. 
UNDERWATER  SEARCH  EXPERIENCE 


several  buses  linking  lower  level  tasks  where  they,  in  turn 
could  also  be  linked  by  buses  to  other  hardware  or 
subfunctions,  etc.  This  hierarchial  design  also  facilitates 
the  implementation  of  Knowledge  Based  Systems  (KBS)  .  Choosing 
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the  optimum  buses  becomes  extremely  important  to  overall 
system  integrity  and  speed  of  response. 

Next,  the  specific  type  of  processor  required  can  be 
chosen  to  match  the  particular  levels  requirements  and 
associated  algorithms.  Whether  one  uses  parallel  processing, 
ASIC  chips,  DSP  chips,  transputers,  or  some  hybrid  combination 
would  depend  on  the  cost,  speed  and  expandability  required  at 
that  level.  In  general,  a  special  function  chip  could  be  used 
at  the  lowest  level  (dedicated,  application  specific,  stand¬ 
alone,  fast,  etc.)  to  evoke  maximum  performance  while  a  state- 
of-the-art  processor  (80X86,  RISC,  transputer  array,  etc.) 
would  handle  the  higher  level  computations  which  would 
necessarily  be  constantly  changing  with  the  environment  and 
mission  objectives. 

Finally,  the  operating  system  must  be  chosen  to  fully 
utilize  the  hardware's  performance  and  integrate  the  various 
levels  of  software.  Real-time  response,  while  supporting  a 
multi-tasking  environment,  must  be  achieved.  For  future 
growth  a  symbolic  language  must  also  be  supported.  The 
ability  to  support  limited  missions  has  already  been 
demonstrated,  but  the  system  level  integration  and  cooperation 
needed  to  support  a  high  level  of  Artificial  Intelligence  (an 
expert  system)  in  an  AUV  is  still  being  researched. 
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VI .  CONCLUSIONS 


A.  RESULTS 

Variable  structure  control  is  a  highly  viable  control 
methodology  for  an  AUV.  Combining  an  adaptive  portion  of  the 
controller  with  the  VSC  method  yielded  satisfactory 
performance  considering  the  nonlinear  nature  of  the  problem. 

Robust  observers  also  have  merit  in  systems  with 
uncertainties.  In  our  application,  this  technique  has  worked 
reasonably  well  for  a  range  of  speed  and  model  parameter 
changes. 

The  ideal  AUV  architecture  would  most  likely  consist  of: 

•  DSP  or  ASIC  chips  at  the  lowest  level  tied  together  with 
a  local  bus, 

•  special  purpose  processors  implementing  the  higher  levels 
and  tied  together  with  a  smart  bus  (such  as  FUTUREBUS+) , 
and 

•  a  special  purpose  optimized  real-time  operating  system  to 
implement  a  multi-tasked  hierarchical  knowledge-based 
event-driven  control  system. 


B.  FUTURE  RESEARCH 

Several  areas  of  AUV  control  are  open  to  investigation 
and  optimization.  The  use  of  DSP  chips  to  pre-process  inputs 
(sensor  fusion)  is  an  important  area.  The  possibility  of 
using  Neural  Networks  (and  associated  parallel  processing)  to 
handle  signal  processing  and  control  problems  is  a  very 
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current  and  far-reaching  issue.  The  appropriate  communication 
links  (bus  type,  protocol,  network  topology,  etc.)  between 
processors  and  the  various  control  levels  is  a  critical  design 
choice  affecting  the  total  system  performance.  Eventually, 
a  fault-tolerant  system  such  as  a  Structurally  Adaptive  System 
(SAS)  must  be  looked  into  to  satisfy  the  reliability  issue. 
A  Real-Time  Knowledge  Based  System  ( RKBS )  should  be  researched 
as  to  the  implementation  of  artificial  intelligence  into  an 
AUV  control  system.  Finally,  a  MIMO  (Multiple-Input-Multiple- 
Output)  VSC  controller  needs  to  be  designed  that  takes  into 
account  all  of  the  inputs  (depth,  speed,  etc.),  the  outputs 
(actuators)  and  as  accurate  a  model  as  possible  in  order  to 
maximize  the  performance. 
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